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Abstract 

We have constructed an e-finite basis of master integrals for all new 
types of one-scale tadpoles which appear in the calculation of the 
four-loop QCD corrections to the electroweak p-parameter. Using 
^ ■ transformation rules from the e-finite basis to the standard "minimal- 

number-of-lines" basis, we obtain as a by-product analytical expres- 
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sions for few leading terms of the e-expansion of all members of the 
standard basis. The new master integrals have been computed with 
the help of the Fade method and by use of difference equations inde- 
pendently. 



1 Introduction 



In many multi-loop calculations the standard method to calculate physical observ- 
ables is to use the traditional integration-by-parts (IBP) method in combination 
with Laporta's algorithm [1,2] in order to reduce all appearing integrals to a 
small set of master integrals. Once this reduction is completed one is left with 
the calculation of the master integrals, addressed in this work. 
One physical quantity which has recently been evaluated applying these methods 
is the contribution from top- and bottom-quarks to the p-parameter at four-loop 
order in perturbative QCD [3] . A subsequent independent evaluation of the same 
quantity has been performed in ref. [4]. This completes the partial result from 
the so-called singlet term already determined in ref. [5] . For the p-parameter, the 
result can be expressed in terms of 63 four-loop tadpole master integrals. A sub- 
set of 13 master integrals was already required in earlier four-loop calculations, 
like the determination of the matching condition for the strong coupling constant 
as at a heavy quark threshold in the modified minimal subtraction scheme [6, 7] 
and the evaluation of the two lowest terms of the Taylor expansion of the vacuum 
polarization function [8, 9]. This subset of master integrals has been evaluated in 
ref. [10] with high precision using the method of difference equations [1, 2, 11-13]. 
All results relevant for the four-loop calculation of ref. [10] have been confirmed 
with a completely different method (see below) in ref. [14]. Some of these master 
integrals have also been found in refs. [5, 13, 15-21]. 

A different method, based on the idea of the e-finite basis, has been suggested in 
ref. [14]. In this approach one avoids so-called spurious poles in the coefficient 
functions, which multiply the master integrals. These poles may arise in general 
while solving the linear system of IBP equations, when a division by e = (4 — li) /2 
occurs. Master integrals which have a spurious pole as coefficient need to be eval- 
uated deeper in the e-expansion. The calculation of each additional order in this 
expansion of a master integral is in general increasingly tedious. In contrast, the 
approach of the e-finite basis exploits the freedom in the choice of master integrals 
in order to select a basis of master integrals in such a way that the coefficients, 
being functions in the space time dimension rf, are finite in the limit e ^ 0. In 
ref. [14] the e-finite basis, shown in fig. 1, has been constructed for the subset of 
13 master integrals discussed above. 

All of them have been calculated semi-analytically. For the problem of the 
four- loop QCD corrections to the p-parameter an e-finite basis with more elements 
is required which is constructed in the following. The corresponding master in- 
tegrals are particularly suited for the evaluation by a semi-numerical method 
based on Pade approximations [14,22,23]. Most of the results are checked in- 
dependently by the method of difference equations and have been used for the 
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Figure 1: Diagrams of the subset of 13 master integrals of the e-finite basis 
addressed in ref. [14]. Note that the last four diagrams are known completely 
analytically and are not considered in the construction of the e-finite basis. 

evaluation of the p-parameter in ref. [3]. They are in full agreement with those 
recently published in [4]. 

This paper is structured as follows. In the next section we present the results 
for the master integrals in the e-finite basis, calculated by the Fade method. 
Then, in section 3, we give results for the master integrals in the standard basis 
evaluated by difference equations. In this case we also include the analytical 
information obtained in section 2. Our summary and conclusions are given in 
section 4. 

2 The master integrals in the e-finite basis 

For the construction of the e-finite basis we follow the algorithm described in 
reference [14] . In the evaluation of the p-parameter we have again a vast choice 
in defining the e-finite basis with elements suitable for the Fade method. We 
exclude 14 diagrams from the construction of the e-finite basis, which are known 
to high orders in e analytically or even completely analytically. Four of them are 
depicted in fig. 1 (T4 i, T5 1, T5 3, Tg 3), the remaining ones are listed in appendix A. 

The application of the method [14, 22, 23] to the actual calculation of the 
e-finite integrals is straightforward. The only adjustment necessary arises from 
the appearance of massless cuts inherent in many of the relevant diagrams. In 
addition to the function used to subtract the high energy logarithms, a second 
function was introduced to subtract the logarithms appearing in the calculation 
of diagrams with a massless cut in the low energy limit. This approach was also 
verified by recalculating the diagrams presented in fig. 1 (e.g. T/4); this time 
performing the cut in such a way that a self energy diagram with a massless 
cut arises. Full agreement between the two approaches was observed. Following 
this procedure we again calculate all divergent parts of the diagrams analytically. 
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and the constant part in the hmit e ^ numerically. All the digits of our 
numerical results in eqs. (3)-(20) are valid digits. The analytical information 
of the expansion then allows to study the cancellation of all divergences in the 
physical problem during the renormalization procedure analytically. Just as in 
the problem of the diagonal current, some diagrams are simultaneously members 
of the standard and the e-finite basis. These integrals are discussed in section 3. 
Now we give the results for the integrals which are different in the two bases, i.e. 
the replaced integrals due to the existence of a spurious pole. As in ref. [14] the 
name of the e-finite master integral is deduced by using the name of the original 
master integral, which has been replaced, and the additional superscript letter 
"/" . The loop integrals are defined as follows. 
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with denominators Dj = {pj — m^). The momenta pj are defined by the dia- 
grams shown below as a linear combination of the four loop momenta fcj. For brief- 
ness we set the single mass scale m of the integrals in the following m = 1, thus 
rrij G {0, 1}. Dashed(solid) lines in all the diagrams denote massless (massive) 
propagators. The integration measure is given by [dkj] = e"^^^ / {i tt*^/^) d'^kj, where 
7£; = 0.577216 ... is the Euler's constant. In this conventions the one-loop tad- 
pole reads: 



Ti = -e^^^r(l-rf/2) 
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The final expressions for the remaining master integrals in the e-finite basis read: 




3C3 1/91 4 15 4 4 

= — ^ TV G log (2) 

2e2 e V360 2 3 ^ ^ ^ 

+ - vr^ log2(2) - 32 Li4 (i) ) + 43.4562683 + 0{e) , (3) 
3 / 
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23.2862341 + 0(e) 



(4) 
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24.8172810 + C(e) 
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- 18.9284641 + 0{e) 
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-5.1732238 + C(e) 



2e 



(5 Cs- 3 Ca) - 1-589373 + 0(e) 
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-8.816973 + 0(e) 
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+ 6.2801043 + C(e) 
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+ 5.68333946 + C(e) 



- 4.872849 + C(e) 



(19) 



(20) 



with the Riemann zeta function (n and the poly- logarithm function Lin{z) being 
defined by: 

^n = Y.i-^ Li.(^) = E7:^- (21) 



k=l 



k=l 



3 The master integrals in the standard basis 

The results from the previous section can also be translated into analytical results 
for the standard "minimal-number-of-lines" basis for a few leading terms of the 
e-expansion by using the IBP-relations between both bases. Since the order of 
the e-finite master integrals in section 2 is known numerically, one can also derive 
one additional order numerically through the IBP-relations for the members of 
the standard basis from the e-finite one. On the other hand one can perform a 
formal expansion of the standard master integrals: 

oo 

T,,= J2 ^T^^ (22) 

^ — ^min 

and insert it into the IBP-relations, which express the e-finite master integrals 
of section 2 in terms of the standard ones. Here we denote each of the master 
integrals according to the following rule: after a capital letter "T" we write the 
number of lines (index i) in the given diagram. The second number (j) enumer- 
ates the different topologies with the same number of lines. If a particular order 
k of the e-expansion of the master integral Tj j is considered, it is denoted by an 
additional last index [k). After performing the e-expansion of the IBP-relations 
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one compares then the different orders in e and obtains thus a hnear system of 
equations. Its solution gives additional relations among particular orders of dif- 
ferent master integrals. They are given in appendix B. 



Another very convenient approach for evaluating master integrals with eight 
or less number of lines is given by the numerical solution of difference equations [1, 
2, 11-13]. This method has been applied directly to the diagrams of the standard 
basis providing an independent check of the results in the previous section and 
leading to much higher numerical accuracy. The idea of this method is described 
in detail in ref . [1] , we restrict ourself on a brief sketch. In this approach one of 
the massive propagators is raised to a symbolic power x 

J{x) = [ [dh][dk2][dh][dh] ^ , (23) 

J JJ1JJ2 . . . JJNd 

and the IBP-method is used to construct difference equations for the master 
integrals by solving a linear system of equations. The appearance of two variables, 
the space-time dimension d and the power x, complicates the solution of the 
system of IBP-idcntities. It is solved with the help of a modified version of 
the setup described in ref. [24] using Form[25-27] and Fermat[28]. In the case 
of the p-parameter also master integrals with increased powers of a propagator 
appear. Then the increased power is chosen in such a way that it is carried by 
the propagator Di, which carries in the above approach the symbolical power x. 
In this case the integral is obtained as a by-product of the solution of the integral 
without extra power. The difference equation for J{x) can be written in the form 

^p,(x)J(x-j) =EE^f , (24) 

j=0 n=l j=0 

where pj and g]"^ are polynomials in x and d. K^'^\x) are sub-topologies of J(x), 
i. e. topologies with at least one propagator less, which also satisfy difference 
equations. At this point the solutions for the sub-topologies are assumed to be 
already known. 

The solution is split into several parts determined by using a factorial series 

as an ansatz. Prom the difference equation one gets the values of /ij and Ki and 
a recursion formula for the coefficients a!i\ This leads to a particular solution 
Jp{x) of the inhomogeneous equation and to one or more homogeneous solutions 
J^\x). The full solution is given by 

J{x)^J,{x) + Y.^,{d)J^^\x) (26) 



6 



with functions rik{d), fixed by boundary conditions. The rik{d) can be determined 
from the low energy expansion of the self-energy diagram which one obtains 
by cutting the line that carries the power x in the diagram J{x). To obtain a 
numerical result, the coefficients a^s^ are calculated as an expansion in e = {A—d)/2 
up to a pre-defined depth emax- Summing the series to a specified s = Smax and 
estimating the remainder gives the numerical result. 

This method provides a possibility to calculate the e-expansions of master 
integrals with high numerical precision up to high orders of e. 

Integrals that we did not calculate by difference equations have been con- 
structed from the e-finite basis. The lowest order, which is only known numer- 
ically, can also be obtained through the IBP-relations from the e-finite basis, 
calculated by the Pade approach, but with less precision. The precision obtained 
with the Fade approach is marked in the results with an underlined (_) digit. For 
11 integrals (T7,io, Ty^n, T7,i2, ^7,13, 7^8,3, Tg^s, Tg^a, T9,4, Tg^g, Tg^e, ^9,7) more precise 
numerical results have in the meantime been obtained in ref. [4]. They are ap- 
pended in the following to our results and separated by a vertical bar (|). All 
given digits of our numerical results are valid digits. The results for the master 
integrals read: 
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+ 80.8955061678534024741104898217973159065185099065309 e 

+ 1085.28365870727983857702515746382256304748449190947 

+ 4545.30388413442580391360145657533885328523257544774 

+ 35998.9938326326556313329606141017307302536969621971 

+ 134897.307552053604704118477425872393007397735113520 + C( 
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161 43 2_,C3 3 /- 

H 71 H V0S2 

192 144 3 2 

+ 7.38924322701052496152773911803203314613623175263798 e 
+ 101.9115108195622938461857471165772655659156888403276^ 
+ 413.971474923722132422941214144351086349919086346111 
+ 2899.99385888238845225724197005268821906921264538028 
+ 6906.33326400669495615717643151117125986689939048417e^ 
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1 19 1 /103 7r2\ 1 /341 19 . r- 4 \ 
= ? + 47 + ^5 (— + T j + ; (ifT + 12 + ^ - 3 «V 
+ 56.2803016207811654954905101988728760185938023940072 

+ 0.62861591651881519024134469355376197437836304764183 e 

- 713.454611137588809603443079328921106975591183151813 

- 3252.49900904808028994099278455214180709190379330122 

- 14732.7748604549077668424518963805693923717102434909 

- 48411.3691880431214632694555990972265902582845398420 + O(e^) , 

(29) 

1 25 1 /257 1 /^433 25 o . /tt 10 , \ 

+ 13.5916504314381310784889154404088589280071858101489 

- 595.0727652159571293998696605616538972651149790881256 

- 6875.84113744124081014720730186766512655050693091046 

- 31308.1849304592520703889751908067773878059674361777 

- 192167.397958807277872297783881991122229270292735722 

- 729545.659023744499039120484965719141090831740292666 + O(e^) , 

(30) 

= l27J + 475 + l27j(^3+rj 

e V 4 12 18 ^7 

- 37.4574649312968260550782554633103343746703536006083 

- 296.071923695141317174869536814657610506616897735608 e 

- 2022.20010652903000709521858304176327369282045103655 

- 8244.01633827202034423181172211035421212582682534110 

- 46235.4396348606702002905412792127287096490637772429 

- 162640.9428015140761391308436008496190801752496553876^ + O(e^) , 

(31) 

7 43 1 / ^ 7 2\ 
= 12^+12^ + 4^ r + 9^ j 
1 /187 43 2 r 7 \ 

+ - + — 77^ + 6 V3S2 + —C3 

e V 12 36 18 ^7 

- 51.2252212801189152418819912274965922080241128506972 
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- 477.947460995194619277285562191897453203074025533346 e 

- 6127.93433249128299350151222877004479343657194425540 

- 23752.7464620411276168084881606893540782888293578810 

- 161921.400195144895866566378264214720989222786416151 

- 555837.682068927549027491 134693752433283969601726090 + 0{e^) , 

(32) 



- 57.5851192392905729834630220283359182548958203213894 

- 226.906688638637087111211732760488182253117551056418e 

- 2357.42884823258106273364123808844182706228997262473 

- 6885.09836097771127634712790473221183461334003969682 

- 51888.4446361973078707935677603480403504827064470340 

- 139706.649615179144594349102097042435381107174124316e^ + O(e^) , 

(33) 



- 59.1987779967727811953138849540094767971219526012487 

- 549.042726795028324060762834253615726678116981354786 e 

- 4293.52617296657045021413760744102602226521264919104 

- 20846.8297417686461476497501335647508755198702169076 

- 110899.174077030571795400409124084216248915113532175 

- 470220.656770151179792855460983971995962315522021662 + O(e^) , 

(34) 



- 62.9263953267520790182858402690009143386475645244229 

- 258.692941139129800870656513171840779122656511379367e 

- 2245.88328900424600912128374855476088078582690149658 
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- 7643.28441760863244533368790482441514380733014950851 

- 48359.2665176408023251619751028356251676249307751892 

- 156259.773282065061616264607990167056276913844653462e^ + 0(e^ 

(35) 

_ 1 5 1 / 7r2\ 1 /67 _ 37 

~ 12^ ^ 12^ ^ 12^2 V ^TyelVY^y^y^^ 

235 65 2 7 4 293 , 

vr^ TT^ H Cs 

4 36 40 18 ^ 

e / 1497 601 , 511 . 1555 , 181 n , 1813 , 

+ -K tt'* H Cs H 71" Cs H Cs 

2 V 2 18 180 9 ^ 27 ^ 15 ^ 

- 3088.55385540997021022612585452885699551634307800797 

- 14498.1092918745193546816622992228822221794484130674 

- 75204.1336314891447627178739121420142533755186792735 

- 315844.168768264796655463071293944512650159761067379 + O(e^) 

(36) 
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^+T7^ + ^ + J+7 ^-6r+ 4^" "3^V 
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- 65696.8989733875039196986140933337667282094540485106 

- 430175.281094824821441568751585189206237721930681617e^ + O(e^) 

(37) 

5 33 1 /775 _ . 

1 /1389 23 2 r 7 . 

+ - + — TT^ + 3 V3 S2 - - C3 

e V 32 8 3 ^ 

+ 147.807870353419457006020027881277752020796006112323 
+ 467.918626785158952385003318864377551570797031294221 e 

- 1988.53979380180044626102894892460888176325158879271 

- 4066.05004959291604399936393129691406383824993650381 

- 84210.5386674689601672799852884017698624685352637357e^ 

- 187546.771892814733657835744490067583363892471336553 + 0(e^ 
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2e4 16 e3 e2 ^ 32 3 J e \ QA 16 
+ 83.1358882000807808465322734822592975462100986465128 

+ 27.1740919045453848408661874005706354672639669874261 e 

- 1554.56495046662560630454485042040757181079394589806 

- 9710.65466554205477699849510014536535786834553811646 

- 59651. 1029754541134802977172410309014967516228078574e^ 

- 263782.712760320335062339831118979031401524684903280 + 0(e^ 
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1 /977 121 2 17 \ 1831 661 tt^ 11 tt^ 187, 



3e V 64 48 3 ^7 128 288 120 1^ 

/42485 335 o 121 , 1327 , 77 2 . 172 , 

- e vr^ H TT^ Ca ti" Cs Cs 

V 256 576 240 36 ^ 27 ^ 3 ^\ 

- 1308.40045975490723588094289794752583132729289225400 

- 6297.44410007372169525480920361163851686356380177551 

- 39381. 3291846386139130702622754006165033230617765764 

- 156879.079031917198094441108322254183097437076390786e^ + C»(e^ 
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1 /175 
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4 60 V 3 --^y 



e V 4 

+ 188.248493841520722512481676436388871071949826368957 
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+ 2625.62618545640179663793065740741631915228108024570 
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+ 280703.977577876240727702494488542000763547356431605 + 0{e^) , 

(42) 



+ 1040.51973958438238064321813618885698154945163332335 e 

+ 2220.52649830986620272733311300227293842057928920038 

+ 19009.3819367933781054485906380197274398245100237012 

+ 29872.2890554813361688482086053842835376892923973047 

+ 325486.618525642205718848428958393877946064518753281 + C(e^) , 
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~ 6e4 + 2e3 + e2 V 3 ^ 18^ 2j 



+ 192.143949985189792783133357258987262614727056122521 

+ 895.583232707441960151372330063515002337522803534682 e 

+ 3015.50609158685899428979534043459900694651565328232 

+ 15431.7141565054206165629974866237765016909275506281 

+ 45950.3851577110455999824698200183839583728323732333 

+ 256376.553671657904111860999953130364103694997971799 + O(e^) , 
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+ 210.851234758305984612090112605434985732961189359568 
+ 790.910352868677999468381955470293781862809895707932 e 
+ 3471.90742617992105936864664247656869003732285384889 
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+ 13388.2008481237206516463529137761462371387930576955 

+ 54666.05398501267487203612268033161919935875865413776^ 

+ 219636.910846771617661122733716410439039757597225140 + C(e^ 
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+ 1861.12230255|2258601194625502 

I - 8297.001296673572071195678072 + C(e^) , (47) 



2 . 1 /13 . 7 _^ , ^. 

'3 



12e^ + 3l^ + ? IT + 36"' + ^^ 



1 /71 14 2 TT^ 68 , 

H ^ TT^H ^ Cs 

3e V 2 3 20 3 ^'^ 

+ 154.604631|544617018499470855434 
I + 894.807380418881266203622406485 e 

I + 2607.13652650580177242514674361 + C(e^) , (4J 



1 1 7 / 1 2\ 1 /7 7 2 5 

= 7 H ^ H ^ l + -7rM TT^ 

12 3e3 12 V 3 y 3e V2 3 3 

+ 10.00336626|15671539152667896043 
I + 121.389726250269120381127113931 e 

I - 338.169587547191526007167931728 + 0{e^) , (49) 
1 3 If 7r2\ 1 / 35 , 

13 



1039 




55 37 4 94 



-C3 + 2V3s2 + 10s2 



+ 



6 



+ 736. 116832299866041378628875269512201227568815410747 e 

+ 3742.45984331167422662482544777380447491847387617756 

+ 11749.4302169340506083438761905399834852011067964612 

+ 62286.5836001930375873157912314697672224972850370815 

+ 184394.067587123000354685027809174508736973039378933 + 0{e^) , 

(50) 



+ 895.267223034507792336545934746971903090074177178563 e 

+ 2851.24487745129920226056389168977151286787950149531 

+ 15769.2924893036658768261052813143569858442259460015 

+ 44094.4294409565473747922022723147055746208297110835 

+ 262762.378299443371029842954614853325970041327609194 + C(e^) , 

(51) 



+ 30.0908437821524694782206570213720645826480229016466 

- 36.8230691904535254355540493262162846080159401689486 e 
+ 260.507827461469390827561610096003346014392421090678 

- 562.809646826601965896917933562966364880244022357921 
+ 2222.67588565840309456690342209333550093550969179227e^ 

- 6072.20564531822489618544102278797200412026771890767e^ + O(e^) , 

(52) 



- 20.8585976739540823192137169923590025170233079442547 
+ 148.814327289184790880960752283863066723514605901922 e 

- 657.418920712168196526528768117001266054030349661525 
+ 3182.47221782449086620274231595579830546198641810013 

- 13795.8211363763697941825608377160301474614616561880e^ + O(e^) , 

(53) 



+ 
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e 

- 19.11186570594|90732311103013898 
+ 141.144175505|094858012998872298 e 

- 605.43119057|5217055408942299967e2 

I + 2945.58041293149149762720437022 + C(e^) , (54) 
e 

- 19.0793750927079960314257405875623941525789104197843 
+ 141.252248186107747164092632797489850337070843975841 e 

- 605.02962120187025890466267145520655064567103503240 
+ 2946.48740435117339409606564510027475583646652772 

- 12664.231982725689641172668688917148051149093738 

+ 54825.71931364361072956312742419885158634286091 + 0{e^) , 

(55) 

= -Cs - 18.76997269185410818615637417243 

e 

I + 142.191941064595949485071666672 e 
I - 600.255511127362803430335175869 

I + 2957.185588135948680830606207646^ + C(e^) , (56) 
e 

- 7.85876279518922242360822990783004542937882100674500 
+ 128.531386676431303444223252316571121713818949997681 e 

- 349.879176215675436895102723350659482216370181950816 
+ 2332.83399433790756728143973285449385093610438970117 

- 8346.39781820188664641354322284607837497614943630437e^ 

+ 40158.5604096624536390379365615442757755373577157863 + C(e^ 

(57) 

e 

- 13.7206302625920162070487342658557388782474967767706 
+ 135.928906108952872038638377087143640185611677958787 e 
-497. 764644552772376477970102393026516204883127610554 



15 



+ 2695.35125395194610632909729485139956434354033178137e^ 

- 11004.9478850481681338765147169019923471294985485170 

+ 49153.9490476063842090788328486862743934773519776031 + O(e^) , 

(58) 

e 

+ 0.374844191926851274910805165871097566997942271321280 
+ 141.683133328263640451220557622143025514756763189200 e 

- 146.6840147856071 1271 1629455842851698248300684175120 
+ 2448.77867872444787260224256254958000162624856674133 

- 5421.70587793103981619383560960628718332878155020120 

+ 38306.6106494119866058980990395532741466206014427522 + C(e^) , 

(59) 

e 

+ 7.30445968508659934705130294982328536505673135057504 
+ 141.443417856878944115040221346329269896865316839719 e 
+ 40.4773793831579631350040651321790566085488080065966 
+ 2244.39062716662154526526360631628245537948712203304 

- 2188.78999525984793680680710077425046912668999350376 

+ 31134.4838067073617857861615173845872142049443337933 + C(e^) , 

(60) 

3 , 1 /tt^ 15 ^ 2^ 

+ 40.4773901150829247806609147361914210374485307884608 

- 100.628977840966481693832556085825783093037821912010 e 
+ 679.056339067891403282049708893797069643756722410073 

- 2387.60452082758635632061839710152233583503731592210 
+ 11259.2595068584556862116342453958677140933294272601 

- 44413.9549636399070618153357715272931135284827413638 + 0{e^) , 

(61) 

- 2.4269563|9537700735 
+ 4.01554669961524192 e 

- 43.6962533603324647 
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Tg 4 




I + 128.936157875347890 + O(e^) , (62) 

+ 0.47361 11472272364450 
I + 1.09585342206826990 e 
I + 5.37764333252884269 

I + 8.82896457590640998 + O(e^) , (63) 

- 3.71140264|536682392682628965373 

I - 2.11520599545877264756135261804 e 
I - 71.9899451389829491830674629550 

I + 41.1881294174309244417864419468 + C(e^) , (64) 

- 6.7284705600856|8105547188977521 

- 26.08764659996|66155389659770717e 

- 214.64771 79124] 11362028052727052 

- 613.715203096|626075654874908838 + O(e^) , (65) 

- 3.4497511|317390349922288 

I + 6.3127694885459812824115 e 
I - 63.668771344187502234181 



I + 196.34402612627359322923 + O(e^) , (66) 

with the symbol S2 defined by: 

S2 = Ch = Im Li2(e^i) = 1.0149416064096536250... . (67) 

In ref. [4] the set of master integrals has been chosen slightly different. Instead 
of the topologies T^^g and Ty^n the topologies Ts^n and Tj^n have been chosen as 
master integrals. These different master integrals can be related to each other by 
the IBP-relations. For these integrals we also give the analytical results, which 
we obtained for the leading e-orders: 




Ts.ll 



4e4 8e3 2e2 V 6 y e V 8 8 ^ , 

+ 30.0253835218317623192659472127405301377841888425075 ' 

+ 254.111031964704093539375826206525134072256004057109 e 

+ 1805.36593366425290012409359932346164787196465139299 
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+ 8898.24672362282959032538518672895187402701187699568 
+ 43751.25518180014447176251 15350236214456071856425684 
+ 193367.012317936435653691364175466657628564760490824 + 0(e^ 

(68) 

1 ^45 7 _2 , 4 , 

b3 



1 


3 






8e4 


" 4^ " 


y 


301 


27 2 

H vr - 

8 


720 


55 


~8~ 





4 12 3 



iS2 + 7^2 



- 17.78744259404|76975226672333232 e 

- 175.317651447|398940503032885592 
+ 313.98310535|3518264983361554717e3 

I - 2142.79466723532874375769322841 + O(e^) . (69) 

All results are in full agreement with those published in [4]. 
An overview over the different members of the two sets of master integrals dis- 
cussed in section 2 and 3 is given in table 1. 

Finally, we want to make a few comments about the results for the master inte- 
grals in the standard basis as collected in eqs. (27)- (69). 

• For the calculation of the p-parameter the coefficients of the master inte- 
grals shown in eqs. (27)-(69) suffer from not more than a 1/e^-pole. On 
general grounds it is clear that in any foreseeable four-loop calculation one 
would always need the analytical terms, displayed in eqs. (27)- (69) of the 
e-expansion of a master integral and additionally up to at most the order 

for the p-parameter. Thus all other higher order e-terms, shown in eqs. 
(27)-(69), will become physically relevant only for five and more loops. 

• The use of the e-finite basis allows for an easy and automatically analytical 
determination of a "trivial" part (that is eventually expressible through 
the three-loop tadpoles) of the e-expansion of a four-loop tadpole master 
integral[14]. This is just enough to control the cancellation of the UV poles 
during the renormalization procedure in completely algebraic way. In order 
to have the final physical result completely analytically one needs to know 
at most one more term in the e-expansion of the members of the e-finite 
basis in an analytical form. In relatively simple cases this could be achieved 
with analytical summation of auxiliary series constructed with the help of 
either the differential relations [29-32] or the difference ones [1,2,11-13]. 
For recent progress in this direction see e.g. [18,19,33,34] and references 
therein. 
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bet oi master 


71 /T 

Megl = 


- i-fe,!!! 


-^7,5, 


-' 7,7, 


rp 

-'7,8, 


rp 

-'7,9, 


-'7,10, 


rp 

-'7,12, 


integrals, which are 






rp 

-' 7,16, 


rp 

-'8,1, 


-'8,2, 


-'8,3, 


rp 

-'8,4, 


-'8,5, 


pnn 1 in nnf Vi V~t/^ g;i c 

v^^JLIOjI 111 UKJkjLI UOjOIO 




T 

8,6i 


T 

^ 8,7, 


8,8, 


T 

^ 8,9, 


T 
8,10, 


T 
9,1, 


9,3, 






^9,5) 


^9,6, 


7^9,7} 










bet OI additional 




" {-'5,2, 


-'5,4, 


rp 

-'5,8, 


-'5,9, 


-'5,10, 


rp 

-'6,1, 


rp 

-'6,2, 


members ol the 




-'6,4, 


-'6,9, 


-'6,10, 


rp 

-'6,12, 


rp 

-'6,13, 


-'6,14, 


-'6,15, 


Git a n n i^rn cji c 

O LOjIILIcIjI Ll UOjOIO 




6,16, 


T 

6,17, 


T 

6,18, 


T 

-'- 6,19, 


-'7,1! 


7,2, 


-'- 7,6, 






^7,11, 


^7,14, 


^7,15, 


^9,4} 








Set of additional 


71 yf 


1^5,2! 


-'5,4! 


-'5,8! 


-'5,9! 


-'5,10! 


-'6,1! 


rpf 

-'6,2! 


rnprntiprs of thp 




6,4? 


6,9! 


6,10! 


6,12! 


-^6,13! 


-^6,14! 


-^6,15! 


e-finite basis 




-'6,16' 


-'6,17! 


-'6,18! 


-'6,19! 


-'7,1! 


-'7,2! 


-'7,6! 






-'7,11' 


-' 7,14! 


-' 7,15! 


^9^4} 








Set of master integrals, 




-- m,l. 


T5,l, 


^5,3! 


^5,5! 


^5,6! 


T5,7! 


^6,3 1 


known analytically to 




^6,5, 


2^6,6, 


^6,7, 


2^6,8, 


77,3, 


2^7,4, 


^9,2} 


sufficient orders in e 



















Table 1: The standard basis is given by Mg^/ U MgtcL U Mana- The e-finite basis 
is given by M^qi U M^pj U Mana- Note that the simple master integrals in the set 
Mana have not been replaced in the case of the existence of a spurious pole in the 
coefficient function, since these master integrals are known analytically to high 
orders in e. The master integrals of the table, which have not been given in this 
work can be found in ref. [10, 14]. 

• One should, however, keep in mind that the numerical accuracy of the 
Laporta method of computing master integrals is huge (30-50 digits!). This 
means that the availability of a completely analytical result is absolutely 
irrelevant for phenomenology. 

4 Conclusions 

We have constructed an e-finite basis for problems related to the p-parameter. 
The calculation of the elements of this basis via the Pade method allows to 
obtain a lot of exact analytical relations among various terms of the Taylor ex- 
pansion in e of the master integrals. In particular we found analytically the 
pole parts of the master integrals in both bases, the standard and the e-finite 
one. For some members of the standard basis analytical results are also ob- 
tained for orders beyond the pole part. With the help of the e-finite basis the 
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cancellation of all divergences during the renormalization procedure can be done 
analytically. These relations are indispensable for any attempts to provide, say, 
the four-loop QCD contributions to the p parameter in a completely analytical 
form. They also serve as a powerful check of available numerical results. We find 
full agreement with the results of the standard basis also calculated in ref. [4]. 
The conversion between the two descriptions as well as all the results presented 
in this work will be made available in computer readable form under the URL 
http : //www-ttp . physik . uni-karlsruhe . de/Progdata/ ttp06/ ttp06-30. 
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A Master integrals excluded in the construction 
of the €-finite basis 

In this appendix we give for completeness the result of the integrals, which have 
been excluded in the construction of an e-finite basis, since they are known ana- 
lytically to higher powers in the e-expansion. The results read: 

; =3e^^^^r(5-2rf) r(3-^) P^d - l) , (70) 



-3e^^^^ (ci-4)P(2-d) P(3-fd) P(3-f) P^(i-2) , (71) 



32e^^7.p(2-d) p3(3-f) P(f-2) 

{d-Af{d-2f ' ^ ' 



8e4^^^P(6-2d) P(6-f) P'(3-f) T^d - l) r(f -5) 
{d-Af{d-2)T{A-d) p2(d-2) ' 

(73) 
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^6,6 



76/ 



7(5,1 



77,3 



T7, 



i2e^^^^r(8-2rf) r(4-d) r(3-f) r(2d-8) r^d-i) 
(rf-3)r(rf-2) r(f -2) 

256 6^^-^^1(6- 1 d) r^d) r3(3- f) r(|d-5) 

(rf-4)3(d-2)3r2(rf- 1) ' 

32e^^T«r(8-2rf) r3(3- f) r(2c/-8) r5(| -l) 
(rf-4)2(c/-2)r3(rf-2) ' 



+ ^ + ^(117 + 7.^)+^ (213 + 28.^ + 10(3) 



, (74) 



(75) 



(76) 



(77) 





1 




12 






+ 


180 




e 


+ 




270 




e2 


+ 




1620 



— (5655 + 1365 + 64 vr^ + 800 Cs) 
180 



- 480195 + 118755 + 22626 vr^ + 1093 vr^ + 273600 Cs 
+ 16800 vr^ C3 - 40740 C| + 387936 (5) + C(e^) , 



(75 



^ + ttL (5 tt^ + 6426 Ci + 5670 Cs) 



25 7r^ + 35771^3 
n'2 + 32130 CI + 22050 C5 + 7350 C5 + 70875 C7 ) + C(e^) • 



e 378 ^ "'^^ 630 

2 I ooncn I '/QKn^S/- , vnc/c; \ 1 /n/'^2^ (79) 



The integrals (70)-(77) are easily found via a repeated application of well-known 
analytical formulas for a one-loop massive tadpole and a one-loop massless prop- 
agator. A less simple diagram (78) can be extracted from [35-37] while the most 
complicated non-planar one (79) from [38,39]. 
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B Relations among particular orders of different 
master integrals of the standard basis 

In this appendix we present the relations among particular orders of different 
master integrals as discussed in section 3: 



(0) _„(o) 11561 16857r^ , Qy/Sg^ , I66C3 .... 

J 6,13 ^6,18- ;^28 288 4 + 2 + 9 ' ^ ' 

T(%-T(?,^^ + lI^ + l^ + 4V3.. + 2.^-^, (81) 



2T<S - SriJ) = I + 51^ + + 12 + 6.1 - H^, (82) 
2 Tfi), - 3 = ^ + 5 + li^ - 9 v^.. + 12 4 - 51^, (83) 



fi 7^(1) , 7^(0) , m(0) 



2097 241^2 35 4 ^ 



^^(1) _ ^(0) , oTn(O) 

" -^5,9 -^5,10 "T -^6,12 



5435 



247 TT^ TT^ /- 2^43C3 

h 9 V 3 So — 24 s, H , 

64 48 2 ^ 6 ' 

.^(0) _ 9 7.(1) ,07^(1) - 

^-'e.lS ^ 6,13 "T 6,14 -^8,10 ~ 

323 71 TT^ 19 TT^ 95 Cs SGTr^Cg 328 



12 36 9 27 



2716223 3533 TT^ 293 tt^ 467 53 



2592 216 180 ^3 

A si 2585 Cs 46 7r2C3 263 C; 



(85) 



(86) 



(87) 



7^(0) , 7.(0) ^(0) _ 
-^5,10 + -^6,17 + -^7,5 — 

55ilZ + 25Ill! + Ill!-9V3.,-15ij + ^-12I^5, (88) 

192 144 180 2 9 9 ' ^ ^ 

c 7^(1) 01 71(1) 071(1) 

"J -^5,8 ^ -'7,14 — 



9 27 5 ' ^ ^ 
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-^5,8 ^-^5,9 ^ 7,10 "T ^ 7,11 ~ 

3471265 4543 229 vr^ 217 S2 



2592 216 45 ^3 

40 si 893 Cs 287r^C3 ^ 2133(5 .q.. 

1 C 9 

-^5,8 -"-^ -^5,9 ^ -^6,16 — 

8415131 _ 39595 TT^ 1039 tt^ _ 349^2 298si 647rMog^(2) 
10368 864 ^ 180 ^ 3 ^ 3 

_ 512 L. (!) - + 1^ + ^, (91) 

T^i + 15 + 2 + 4 T«o + 2 T^W - 6 T^^ + 2 T^Jg = 

957005 110277r2 _ 23 tt^ 53 ss _ 140 si 
2592 ^ 216 30" ^ 3~ 

238 Ca lOlTT^Cs , 974 Cs 



+ -rr^, (92) 



9 27 15 

+ 15 _ 2 + 4 + 2 \\ - 6 T^i, + 12 T^Jg - 4 r«3 

(1) (0) _ 2290913 14459 TT^ _ 79 7r^ _^ 179^2 

+bJ6,i4 + bJ7,8 - 2592 + 216 45 ^ 

212 si 100 C3 ^ 537r^C3 ^ 634 Cs 

3~ + ~9~ + ^^ + ^' ^^^^ 

/I 7^(1) 1 9 r-p(l) r T-t(O) r\ rp(l) rp(l) . q rp{l) 119 rp(0) 

^-'5,8 ^^-'■5,9 "^-^5,10 ^-^5,10 6,11 "T -^6,12 ^6,13 



964973 6613 43 tt^ 



1296 216 72 



184S2 ^ 52^1^1618(3^ 29^ ^^g^^^ ^^.^ 



x/3 3 9 27 

T^s + 15 + + 4 + 2 ri^A - 6 T^}, + 6 7^53 - 4 7;% 

+6 t(^) + 3 t(°) 3 t(°) - ^^^^^^^ + 29521 TT^ _ 313 tt^ 
+6 76,14 + 3^6,17 3^7,7- ^^g^ + 180 

8s2 176si 2551C3^ 35 7r2C3^ 5372C5 
+7f"^ 1^ + ^^ + ^^' ^''^ 
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^ rp{l) qi T-i(l) 07^(0) 9 T-i(l) rp(l) . n _|_ Ofl T''-'^-' Hl^^^^ 

'-'-5,8 '-'-'--'5,9 "-'-'5,10 ^-'S.IO 6,11 ^ ^6,12 "T 6,13 °-^6,13 

+12 T« +2T(°) + 6T('^)- 590641 2627 vr^ 707 tt^ 
508 S2 ^ 376 si ^ 53 tt^ log(2) ^ 4 tt^ log^(2) 4 log^(2) 



3 ' 15 3 
167rMog^(2) _ 16 log^(2) _ 



^-^8^1- 2127 Cs 
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